ON THE MARKOV TRANSITION KERNELS FOR FIRST-PASSAGE PERCOLATION ON 

THE LADDER 



ECKHARD SCHLEMM 

ABSTRACT. We consider the first-passage percolation problem on the random graph with vertex set IN x 
{0, 1}, edges joining vertices at Euclidean distance equal to unity and independent exponential edge weights. 
We provide a central limit theorem for the first-passage times l n between the vertices (0,0) and («,0), thus 
extending earlier results about the almost sure convergence of l n /n as n — > °°. We use generating function 
techniques to compute the n-step transition kernels of a closely related Markov chain which can be used to 
calculate explicitly the asymptotic variance in the central limit theorem. 



1. Introduction 

The subject of first-passage percolation, introduced in Hammersley and Welsh (1965) in 1965, is the 
study of shortest paths in random graphs. Let G = (V,E) be a graph with vertex set V and unoriented edges 
E C V 2 and assume that there is a weight function w : E —> R + , For vertices u, v g V, a path joining u 
and v in G is a sequence of vertices p u ^ v — {u = po,p\,. . . ,p„-i,p» — v} such that (p v ,p v+ \) g E for 
< V < n. The weight w(p u _^. v ) of such a path is defined as the sum of the weights of the comprising 
edges, w(p u ^ v ) := £y=nw((/?v>.Pv+i))' The first-passage time between the vertices u and v is denoted by 
do(u,v) and defined as dc(u,v) := inf{w(p),p a path joining u and v in G}. 

First-passage percolation can be considered a model for the spread of a fluid through a random porous 
medium; it differs from ordinary percolation theory in that it puts special emphasis on the dynamical aspect 
of how long it takes for certain points in the medium to be reached by the fluid. Important applications 
include the spread of infectious diseases (Altmann (1993)) and the analysis of electrical networks (Grim- 
mett and Kesten, 1984). Recently, there has also been an increased interest in first-passage percolation on 
graphs where not only the edge-weights, but the edge-structure itself is random. These models, includ- 
ing the Gilbert and Erdos-Renyi random graph, were investigated in Bhamidi et al. (2010); Sood et al. 
(2005); van der Hofstad et al. (2001) and found to be a useful approximation to the internet as well as 
telecommunication networks. 

Usually, however, the underlying graph is taken to be Z 2 and the edge weights are independent ran- 
dom variables with some common distribution P, see e.g. Smythe and Wierman (1978) and references 
therein. Interesting mathematical questions arising in this context include asymptotic properties of the 
sets B t := {u <E Z 2 : d Z 2(0,u) < ?} (Kesten, 1987; Seppalainen, 1998), and the limiting behaviour of 
c/ Z 2((0,0), (n,0))/n. The latter is known to converge, under weak assumptions on P, to a deterministic 
constant, called the first-passage percolation rate; the computation of this constant has proved to be a 
very difficult problem and has not yet been accomplished even for the simplest choices of P (Graham 
et al., 1995). An exception to this is the case when the underlying graph G is essentially one-dimensional 
(Flaxman et al., 201 1; Renlund, 2010; Schlemm, 2009). 

In this paper we consider the first-passage percolation problem on the ladder G, a particular essentially 
one-dimensional graph, for which the first-passage percolation rate is known (Renlund, 2010; Schlemm, 
2009). We extend the existing results about the almost sure convergence of dc((0,0), (n,0))/n as n — > °oby 
providing a central limit theorem as well as giving a complete description of the n-step transition kernels 
of a closely related Markov chain. Our results can be used to explicitly compute the asymptotic variance in 
the central limit theorem. They are also the basis for the quantitative analysis of any other statistic related 
to first-passage percolation in this model. In particular, knowledge of the higher order transition kernels is 
the starting point for the computation of the distribution of the rungs which are part of the shortest path. 
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Figure 1. The ladder graph G„. The edge weights X v , Y v , Z v are independent expo- 
nential random variables. 



The ladder model is worth studying because it is one of the very few situations where a complete explicit 
description of the finite-time behaviour of the first-passage percolation times can be given. 

The structure of the paper is as follows: In Section 2 we describe the model and state our results; Sec- 
tion 3, which contains the proofs, is divided in three subsections. The first is devoted to the central limit 
theorem, the second presents some explicit evaluations of infinite sums which are needed in Section 3.3, 
where the main theorem about the transition kernels is proven. We conclude the paper with a brief discus- 
sion. 

We use the notation 8 P ^ for the Kronecker delta and & p q as well as & p >q for versions of the discretized 
Heaviside step function, precisely 



J p,q 



1 if p = q 
else 







p,q 



1 if p > q 
else 







1 if p < q 
else 



The symbol (k) ] stands for k\{k+ 1)!. We denote by R the real numbers and by Z the integers. A subscript 
+ (— ) indicates restriction to the positive (negative) elements of a set. The symbol y stands for the Euler- 
Mascheroni constant and E denotes expectation. 

2. First-passage percolation on the ladder 

In this paper we further investigate a first-passage percolation model which has been considered be- 
fore in Renlund (2010) and also in Schlemm (2009). We denote by G„ the graph with vertex set V„ = 
{0, . . . , n} x {0, 1 } and edges 

£-„ = {((/,()),(/+ 1,0)) :0 </<«}, 
^„ = {((/,l),(/+l,l)):0 </<»}, 
iT„ = {((;, 0),(/,l)):0</<«}. 

The edge weights are independent exponentially distributed random variables which are labelled in the 
obvious way Xf, F; and Z,-, see also Fig. 1. By time-scaling it is no restriction to assume that the edge 
weights have mean one. We further denote by /„ the length of the shortest path from (0,0) to (n,0) in 
the graph G„, by l' n the length of the shortest path from (0,0) to (n, 1) and by A„ the difference between 
the two, A„ = l' n — 1„ It has been shown in Schlemm (2009) and also, by a different method, in Renlund 
(2010) that \\m n ^ool n /n almost surely exists and is equal to the constant % = | — jj^)' wnere Jv are 
Bessel functions of the first kind. (See Definition 1 or Abramowitz and Stegun (1992) for a comprehensive 
treatment.) This constant is called the first-passage percolation rate for our model. The method employed 
in Schlemm (2009) to obtain this result built on Flaxman et al. (201 1) and consisted in showing that there 
exists an ergodic IR x Revalued Markov chain M = (M„)„>o with stationary distribution ft and a function 



Explicitly, 



R such that 



X = E/(M») = 



xIR 3 



/(m)£(dm). 



M„ 



(A„,X„ +1 ,y„ +1 ,Z„ +1 ) and / : (r,x,y,z) ^ min{r + y + z,x}. 
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Throughout we write m = (r,x,y,z) for some element of the state space R x Rl, In order to better under- 
stand the first-passage percolation problem on the ladder it is important to know the higher-order transition 
kernels : (R x R 3 + ) x (R x R^/j — >• R + of the Markov chain M. They completely determine the dynamics 
of the model and are defined as 

(2.1) K n (m',m)<dm = P(M n e<dm\M =m'), m,m' eRxR 3 + . 

The first result shows that it is sufficient to control the transition kernels K" : R x R — >• R + of the Markov 
chain A = (A„)„>o which are analogously defined as 

(2.2) K n (r',r)dr=P(A,,edr\A = r , ) : r, / G R. 
For convenience we define K°(r' ,r) := 5/(r), the Dirac distribution. 

Proposition 1. For any n > 1, denote by K" the n-step transition kernel ofM defined in Eq. (2.1). Then 

(2.3) K n {m',m)=e- {x+y+ ^K n - 1 (min{/ +/ ,x' + z'} - mm{r' +y' +z',x'},r) . 
Moreover, the stationary distribution ft ofM is given by 

(2.4) n{Am) = e- ix+y+z) <d 3 {x,y,z)n(<dr), 
where 

(2.5) n{&r) = ^-^e-^h(2e-^)&r 

is the stationary distribution of A. 

Next, we state a central limit theorem for first-passage percolation times on the ladder which has been 
implicit in Schlemm (2009) and which was the motivation for the current paper. In Ahlberg (2009); Chat- 
terjee and Dey (2009) a central limit theorem has been obtained for first-passage times on fairly general 
one-dimensional graphs by a different method. The question of how to compute the asymptotic variance 
has, however, not been addressed there. We denote by / the mean corrected function f — %. 

Theorem 1. For any integer n > 0, let l n denote the first-passage time between (0,0) and (n,0) in the 
ladder graph G„. Then there exists a positive constant <J 2 such that 

(2.6) k^*4jv(0,<7 2 ), 

V n 

where N(0, <7 2 ) is a normally distributed random variable with mean zero and variance (J 2 and denotes 
convergence in distribution. Moreover, 

(2.7) C7 2 =/ /(m) 2 7T(dm)+2f; / f{m)P n f{m)fc(Am), 
where 

(2.8) P"f(m) = E [f(M„)\M = ml = / f(m')K n {m,m')dm' . 

JRxR 3 + 

Equation (2.7) shows that in order to evaluate the asymptotic variance of the first-passage times one 
must know the transition kernels K" . In the next theorem we therefore explicitly describe the structure of 
the transition kernel K" and thus, by Proposition 1, of K n . To state the formulas in a compact way we define 
the five functions 

, z-2J 2 (2^) 2 2( Z -1)+2J (2^) 

s {z)= , s (z) = , 

z z 



G(z) 



5-47 + 2,^f^-21og, 



z A 

a{£) = 4 ^J 2 (2 y/z) [V~zh(2V~z) + (z-l)h(2 V~z)} ^ 
z 2 

H(z) -^Y^y [3z - 3 + 2n Y (2 y/z) + J (2 y/z) (5 - 4y-21ogz)] 
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The Bessel functions J v and Y v are defined in Definition 1 and treated comprehensively in Abramowitz 
and Stegun (1992). 

Theorem 2. The transition kernels K", defined in Eq. (2.2), satisfy K n (y ,r) = K n (— r', — r). For r > f/ze 
values K n (r ,r) are gz'vezz fey 

E™=o< ? e prM9+2)r 

(_l)i-l < /-(«+l)r 2 ( _iy vV - P ( r '- r )-« r ^_ ( +2)r , 

(2 - 9 ' S {^Ay- +K=0 (p)\n-p-2)> +Lp.ci=V b P,q e U<r<r ; 

( _ 1)M -I e -(„-l)/-,- ( _i)> Ve - P (/- r )- m - pr /_ (?+2 ), , 

(n-l) ! T ^p=0 (p)!(„-p-2) ! ^^p.^C-p,^ ^ 

where the coefficients a" p q , b" p q and c" p q are determined by their generating functions: 
i) the generating functions A PjI? (z) = Y,n=i a " p qZ n > PA>Q> are given by 

(<z) ! 



(2.10a) A Ul (z) = K -j^r a(z), 



(2.10b) A p>q (z) = 2( & 1 A 1|g (z), p > 2, 

(2.10c) Ao i? (z) = ^A li9 ( z ); 

i — z 

iij f/ze generating functions B /X(? (z) = LJT=i ^p g 2 "' P>9 — 0> flre If've/i fry 
(2.11a) B li? (z) = ^f G(z)-A l!? (z), 

(211b) *™^'^ B '^ + ihThwTTy pi2 ' 

1-z (?)> 
z'z'z'J f/ze generating functions C p , q {z) = L5T=1 c p </ z "> /?i<7 > 0, are given by 



(2.11c) b ^(z) = t -^Bi i9 (z) + ^-H(z); 



(-z)«+ 2 



(2.12a) C , ? (z) = d q z? +2 +B , 9 (z) - 

(2.12b) C 1>? (z) = 



S 1(Z) + ZS2(Z) 



2(1 -Z) 



Ai, 9 (z)--C , ? (z), 



(2.12c) C M (z)= V , , Ci, g (z), P >2 

and the numbers d q are determined by their generating function D(z) — L„=o dqZ q given by 

(2.13) D(z) = - [ y/ih (2 Vz) (27+ logz) - it \fzY\ (2 ^/z") - 1] . 

By the properties of generating functions the coefficients a" } q , b' p q and c p q are determined by the deriva- 
tives of Ap i? , B p<? and C p ^ evaluated at zero. These derivatives are routinely calculated to any order with 
the help of computer algebra systems such as Mathematica. Table 1 exemplifies the theorem by reporting 
the values of the coefficients a n p q , b p q , c p q in the case n = 4. Using these results the expression (2.7) for 
a 2 can be evaluated explicitly in terms of certain integrals of hypergeometric functions; the computations, 
however, are quite involved and the final result rather lengthy, so we decided not to include them here. 

3. Proofs 



3.1. Proofs of Proposition 1 and Theorem 1. In this section we present the proofs of the relation between 
the Markov chains M and A and of the central limit theorem. 
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(A) values of 

o 1 2 



(B) values of bt„ 



q 
l 



2 3 



(C) values of „ 



2 3 



115 

96 

J_[ 

36 
11 
108 

J_ 

72 

1 

" 1440 



_ 17 

36 
_J_ 

36 
J_ 

12 

1 

144 



24 
J_ 
12 
J_ 

72 







1 

" 144 







721 
432 

241 

540 

3_ 

16 

1 

216 



_ 13 
12 
17 
48 

J_ 
36 






J_ 

12 

j_ 

36 















721 

432 
241 
540 
1 

144 
1 

216 





_ 13 
12 

12 

48 
_j_ 

36 






18 
J_ 

36 















Table 1. Coefficients of the four-step transition kernel K , as given in Eq. (2.9) 



Proof of Proposition 1. Since 

a„ =/;,-/„ 

=min{4_ 1 +F„,/„-i +X„+Z„}-min{l' n _ l +Y„+Z„,l n - l +X„} 
= min{A„_ 1 +F„_ 1 ,X„_ 1 +Z„_ 1 }-min{A„_ 1 + T„+i +Z n _ u X n _ l } 

it follows at once that Mo being equal to some m' = (r 1 ,x ! ,y' ,z') implies Ai = minlr 1 +y' ,x' + Z 1 } — min{r' + 
y 1 + z*,x!} and thus the Markov property of A together with the independence of the edge weights implies 
that for any integer n > 1 the conditional probability P (M„ e dm|Mo = m') is given by 

e - x+y+z d 3 (x,y,z)P {A„ G dr|Aj = min{r' +y',x'+z'} - min{/ +/ + z',x'}) . 

The homogeneity of the Markov chain A then implies Eq. (2.9) because for n = 1 we clearly have 

P (Mi G dm\M = m) = e _ ^ +> ' +z) 5 min{r / + yy +z / } _ min{/+ y + ,/ j y } (r)dm. 

Equation (2.4) about the stationary distribution of M is a direct consequence of the fact that the edge weights 
X n+ \, Y n+ i ,Z„ + i are independent of A„ and expression (2.5) was derived in Schlemm (2009, Proposition 
5.5.). ' a 

Next, we prove the central limit theorem and the formula for the asymptotic variance a 2 . 

Proof of Theorem 1. We apply the general result Chen (1999, Theorem 4.3.) for functionals of ergodic 
Markov chains on general state spaces. We first note that 

x2#,, x 2J 1 (2)-3J (2)+ 2 F 3 ({1,1}, {2,2,2};-!) -1 
f{m) n{<Am) = — — < oo. 

It then suffices to prove that the Markov chain M is uniformly ergodic, which is equivalent to showing that 
the Markov chain A is uniformly ergodic. We use the drift criterion Aldous et al. (1997, Theorem B), which 
asserts that if there is a sufficiently strong drift towards the centre of the state space of a Markov chain, it 
is uniformly ergodic. Using the Lyapunov function V(r) = 1 — e~l r l as well as 



(3.1) K(r',r) = 



if r' < r < OV/ > r > 



e \ r ' 2r l else 



for the one step transition kernel of A (Schlemm, 2009, Proposition 5. 1) we obtain 
y/(r) :=P l V{r) =1-E,.<T |A ' 



= 1- / e-\P\K(r,p)dp= l - 



3-2e-^+e- 2 ^ 



It is easy to check that 



\ff(r) < V(r) — r^r , \r\ > 1 and sup \jf(r) < sup y/(r) = - < oo. 

10 |,.|<i r g|R 2 



6 



ECKHARD SCHLEMM 



Since the interval [—1,1] is compact and has positive invariant measure it is a small set (Nummelin and 
Tuominen, 1982, Remark 2.7.) and it follows that A is uniformly ergodic, which completes the proof. □ 

3.2. Summation formulas. In this section we derive some summation formulas which we will use in the 
proofs in Section 3.3. Some of them are well-known, others can be checked with computer algebra systems 
such as Mathematica, a few (formulae 8,9, 11 and 1 2) seem to be new. The sums will be evaluated explicitly 
in terms of Bessel and generalized hypergeometric functions, which we now define. 

Definition 1 (Bessel function). Let X be a real number in R\Z_. The Bessel function of the first kind of 
order X, denoted by Jx, is defined by the series representation 

°° (—l) k rx\ 2k +^ 

° 2) h{x) -^mWTxTV)(2) • 

For any integer V the Bessel function of the second kind of order V, denoted by Y v , is defined as 

(3.3) Y v (x) = hm — . 

k^v sinA7T 

It is well-known that Bessel functions satisfy the recurrence equation J v (*) = 2( ' V ~ 1 ' > J v -i {x) — Jv-2{x) 
(Wolfram Research, Inc., 2010, Formula 03.01.17.0002.01), which we use without further mentioning to 
simplify various expressions. 

Definition 2 (Generalized hypergeometric function). For non-negative integers p < q and complex num- 
bers a = ai,...,a p and b = b\, . . . ,b q , bj$~L-, the generalized hypergeometric function of order (p,q) 
with coefficients a, b, denoted by p ¥ q (a,b; ■), is defined by the series representation 

(3.4) p F q (a,b;x = L TjTVm' 

where (z)k denotes the rising factorial defined by (z)k = T(z + k)/T(z). 

In particular we will encounter the regularized confluent hypergeometric functions oFi, which are de- 
fined by oFi ({}, {b}; — z) = jtooFi ({}, {b}; — z); in the next lemma we relate their derivative with respect 

to b to certain sums involving the harmonic numbers := Y,n=l 1 l n - 

Lemma 1. Denote by oFi the regularized version of the hypergeometric function qFi. It then holds that 
i) for every positive integer V, 

(3.5a) A F l({ } ){fo}; _ z ) 



ii) for every positive integer V, 

^ «>,{<*-= 



k=0 



t_0 * ! 



Proof. For part i), we differentiate the series representation (3.4) term by term. Using the definition of 
the Digamma function F as the logarithmic derivative of the Gamma function F as well as the relation 
F(k) = -y-t-Hfc-!, keZ+ (Wolfram Research, Inc., 2010, Formula 06.14.27.0003.01), we get 



<d 1 



and thus 



dbT{k + b) 
A f i({})W; - z ) 



F(k + V) Y-H 



k+v-l 



Y{k + v) (jfc+v-1)! 



which concludes the proof of the first part of the lemma. Part ii) is shown in a completely analogous way, 
using the relation lim^_ > _ m F(ju)/r(/x) = (— l) m+1 m! for every non-negative integer to, which follows 
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from the fact Wolfram Research, Inc. (2010, Formula 06.05.04.0004.01) that the Gamma function has a 
simple pole at — m with residue (— l) m /ml. a 

Formula 1. 



y (-z) k = 21 2 (2Vz)-z 
£*!(* + 2)! 



2- 



Proof. Immediate from the definition (3.2) of Bessel functions. 
Formula 2. 



h (k+lV 2 



1 



k=1 {k+i)i* z 
Proof. This is also clear from the definition (3.2) of Bessel functions. 



S 2 (z). 



Formula 3. 



s 3 (z)= £ 



(-z)"- ei 



n= q +2 ("-?- 2 ) ■("-?) 

Proof. Shifting the index of summation n by q + 2 we obtain 

(- z y+ 2 



= l-h (2y/z)- V~zh(2V~z)- 



s 3 w = E 



„f «!(n + l)!(n + 2) 2 to 



(-z) 



B+l 



£j W ! ' T< „to "< 2 
=1-J 2 (2 v^)- ^zJi(2V5) 
by the definition (3.2) of Bessel functions. 



\n+2 



l-Z 



(« + l)!(n + 2)! (n + 2)! 2 



Formula 4. 



q-2 



El(CZ) i5>£(*) ! (?-*-2) ! (* + 2) 2 



MS 



Proof. By Fubini's theorem we can interchange the order of summation and then shift the summation index 
q by k + 2 to get 

^> L mk+2) 2L {q y. ■ 
By definition (3.2), the second sum equals Ji (2 \fZz) I v?Z an d so the claim follows with formula 3. □ 



Formula 5. 



q-2 



2(-Cz) 4 



2J 2 (2 7ciJi(2yci 



,=oS*!(*+2)!(«-*-2)' 
Proof. We can interchange the order of summation and shift the index q by k + 2 to obtain 

v2 " 2(-g Z )*+ 2 £ (-&)« 

h k\(k+2)\ L {q y. 

The first factor equals 2£,zh (2 -\/Cz) — (C z ) 2 an d me second factor equals Ji (2 \/£z) / both by 



definition (3.2), and so the claim follows. 



Formula 6. 



s 
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Proof. Clear from the definition (3.2). 
Formula 7. 



1-J (2v^)- v^zJi(2v^). 
Proof. This follows from the decomposition ^ + ^, 2 = — an( ^ definition (3.2). 

Formula 8. 

n i M = f { ~ z)k+2 V 21 + 1 

3z 2 /? 
=4- - z- 2 - *zY 2 (2 + -y- Ji (2 Vz) [4y- 3 + 21ogz] 



Proof. First we note that 
(3.6) 



-Jo (2Vz) 



E 2/ ^' 



5z _ 
1 + y -27z-zlogz 



1 1 



I l + l 



2H, 



5£ + 3 
'2{k+l) : 



where denotes the k th harmonic number. The first contribution to U (z) can therefore be computed as 

1 ~ (5k + 3)(-z) k 1 3 1 ~ (-z) k 5 y (-z) k 
2£}*!(* + 2)!(* + l) z 4 z£ (k)'- 2^k\(k + 2)\ 

4Ji (2 Vz) + 10 Vzh (2 Vz) ~ 3z 3/2 - 4 Vz 



4z 3/2 



The other contribution to U 1 (z) is, with the help of Lemma l,ii), obtained as 



of- ttf H _i 



7zJ 2 (2 V5) -z - 1 - ^oFi ({}, {^}; -z) 



By Wolfram Research, Inc. (2010, Formula 07.18.20.0015.01), 

1 



iLoP, ({},{&}; -z) 



[ttzY 2 (2 Vi) - VzJi (2 Vz) [2 + logz] + Jo (2 ^) [zlogz- 1]] 



j=-i 



and the result follows after combining the last four displayed equations. 
Formula 9. 



9-2 



y I \ \" ' ~ ' f 21 + 1 J 1 

^h^- k - 2 y- k] -i k + 2 y-kA i + i ) 



(2^ 



U'(Cz) 



Proof. Interchanging the order of the first two summations and shifting the index g by k + 2 we find that 

r 1 ir \ = f ( ~ Cz) * +2 f izM! f 2/ + 1 



As before the sum in the middle equals Ji (2\/(z) / V Cz an( J so the claim follows with formula 8. □ 
Formula 10. 



S 4 (z)= E 



c-zr 



n—q-yl 



(n — q — 2) - (n — q — l) 2 



:z 2 2 F 3 ({1,1},{2,2,2};-z). 



Proof. After shifting the index « by g + 2 this follows immediately from the definition (3.4) of the hyper- 
geometric function. □ 
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Formula 11. 

tt2/_\ _ y (~z) k+2 y 2/ + 1 

= |[-5 + 3z + 2ttY (2Vz)-2zaF 3 ({l,l},{2,2,2},-z)+Jo(2^)(5-4y-21ogz)]. 
Proof. The proof is analogous to formula 8. Using Eq. (3.6) the first contribution to U 2 (z) is 

(-z) k 



1 " (5k + 3)(-z) k+2 

ity (k+iy 2 (k+i) 



5-3z_5_^ {-z) k y 

2 ^=o k]2 h( k + i v 2 ( k + 1 ) 2 



2z 



(3.7) 



1 

~2z 



[5-3z-5J (2 vS)-2z 2 F 3 ({l,l},{2,2,2};-z) 



where we used the definitions (3.2) and (3.4). For the remaining part we first use Lemma l,i) and the 
identity Ji (2 yfz) / y/z — J2 (2 y/z) = Jo (2 y/z) to compute 



(— z)* 
= £ (£+l)!(fc + 2)! 



Hi+2+ | 1 ^T27! H ^ 2 



3 1 1 

o+ 

2 z z 

3 1 7 



fo^ H * +1 ^fo^T27I H/:+2 



Jo (2 V?) 



^oFi({},W;-z; 



7=2 



■z^oFi({},W;-z) 



i>=3 



It then follows from the explicit characterization of ^oFi ({}){&};— z)L = _ v , V G Z + , given in Wolfram 
Research, Inc. (2010, Formula 07.18.20.0013.01) that 



Using this we get 

-z) k 



1 



2z 3 /2 



[ yfz [2 - 3z + TT Yq (2 Vz)] - 2Ji (2 Vz) - V^Jo (2 vS) [27+ logz]] 



(-z) k 



(-z) k 



(-z) k 



2 ^(^+l)! 2 * &(* + l)! 2 ' +2 &(*+l)! 2 (* + 2) (*+!)!>(*+ 1) 



k=i k^^^j- k=i 
[ yfz [2 - 3z+ % Y (2 Vz)] - 2Ji (2 Vz) - v^Jo (2 Vz) [27+ logz]; 



(3.8) 



k=\ 
1 

= 2zV2 

"1 1 1 Ji(2Vz) 

2 z + z 3 / 2 



•l- 2 F 3 ({l,l} ! {2,2,2};-z). 



Combining Eqs. (3.7) and (3.8) completes the proof. 



Formula 12. 



T 2 ((:7] = y y (-^ X 21 + 1 



h (2 y/Cz) 



u 2 (Cz) 



(W 3/2 

Proof. Interchanging the order of the first two summations and shifting the index q by k + 1 we find that 



T 2 ((:7] = y (-Zz) k+2 y (-tz)«- 1 * 2Z + 1_ 



By definition (3.2), the middle sum is equal to — Ji ^2 \J~£/zj /(£z) 3//2 and so the result follows readily 
from formula 11. □ 
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Formula 13. 



q=0k=l 



2(-£z)« 



Cz-l+Jo (2 v^z)]Ji (2 v^) 



-k-l) [ (k+l)\ 2 
Proof. The same as before so we omit it. 



Formula 14. 

£ 4 (Cz) 



00 q-l 



(~&) S 



k-l)'(k- 

Proof. Using Fubini's theorem we get 



CzJ. (2^Cz 2F3 ({1,1}, {2,2,2};-^) 



(-&)* 



(-Cz)«- k 



The first factor is equal to 2F3 ({1, 1}, {2,2,2}; — £z) by definition (3.4) and the second factor equals 
- VTzJi (2 VTz) by definition (3.2). □ 

Having evaluated these sums we now turn to relations between them which we will also need and which 
are proved by writing out the expressions and straightforward computations. The first one only involves 
the five functions appearing in the statement of Theorem 2. 

Lemma 2. For almost every complex number z with respect to the Lebesgue measure on the complex plane, 



(3.9) 



S 2 (z) 
2(1 -z) 



G(z)- 



«(*)• 



H(z) 
2 



3z z 



= 0. 



's'M , s 2 (z) | r 

z 1— z z_ 

Lemma 3. For almost every complex number z with respect to the Lebesgue measure on the complex plane, 



(3.10) 



U 1 (z) + S J (z) + 



G(z) = lz 2 -z-l. 



3.3. Proof of Theorem 2. In this section we prove the main result, Theorem 2. First, however, we spend 
a closer look on the coefficients a n p q , b" p q and c" p q which are defined implicitly through the generating 
functions (2.10) to (2.12). 



Lemma 4. For any integers p,q > the following hold: 



(3.11) 



S p ,iS q fi, 



0. 



0. 



"p,q "P,i"9,ui "p,q "1 ~p,q 

Proof. This follows from evaluating the derivatives of the generating functions A M , B p q , C Pt q at zero. 



Next we derive some useful relations between the coefficients a" „ b"„ and c" ,. 



These will be the 

main ingredient in our inductive proof of Theorem 2. The general strategy in proving the equality of two 
sequences (s n ) n > 1 and (.?„)„> 1 will be to compute their generating functions £»>i s n z" and Y. n >i ^nz" and 
to show that they coincide for every z. The validity of this approach follows from the well-known bijection 
between sequences of real numbers and generating functions (see e. g. Wilf (2006) for an introductory 
treatment.) We will constantly be making use of the convolution property of generating functions. By 
this we mean the simple fact that if (s„)„>i is a real sequence with generating function S(z) and (f„)„>i is 
another such sequence with generating function T(z), then the sequence of partial sums (Ly=i f n-v*v)n>i 
has generating function S(z)T(z). We also encounter generating functions of sequences indexed by q 
instead of n. In this case we denote the formal variable by £ instead of z and sums are understood to be 
indexed from zero to infinity. 

Lemma 5. For every integer n > 1, q > the coefficients defined by the generating functions given in 
Theorem 2 satisfy the relation 

q—2 n n 



(3.12) 



q-2 un 

V k,q-k-2 



E 

k=0 



-k.q-k-2 

k + 2 



(-1)>-1) 



n-2 



(-1)" 



7l2 



^{k)\n-k-2)\k + 2f 
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Proof. These equations are true for all n > 1 if and only if the corresponding generating functions coincide. 
Multiplying both sides by z", summing over n and using the recursive definitions of the generating functions 
as well as the fact that, by formulae 1 and 2, 



Y (-Z) k = lew, 

£i*!(* + 2)! 2 1 h 
we find that the claim of the lemma is equivalent to 



tl (^ + 1)! 2 



k=l 



21+1 1 



9 -2 

° = ^ l M(k + 2)\(q-k-2y £{1(1 + 1) ' 2 



(~zf 



0, 



9,2 



■ q—2 



2(-l) 



2 ^ fc!(/fc + 2)! 9 

5 3 y 2(- Z )^-3 

4 £ 1 *!(* + 2)!fa-*-2) 



-Z%-2 

(-;)* 



<?! 2 £- (*) ! (?-*-2) ! (* + 2) 2 



where we have used Lemma 2 to simplify the coefficient of the last sum. To show this equality for all 
non-negative integers q we compare the ^-generating functions and must then show that 

(3.13) 0=T\Cz) + ^[t\Cz) + ((S 1 (Cz)-1)D(Cz)]~T 2 (Cz)+^(Cz) + It. 2 (Cz), 

where closed form expressions for 



9-2 



(-Cz)« 



9=0 £=0 



^{ky(q-k-2)\k + 2r 



9=C 

T 1 ^) :=£ 



^ 2 2(-Cz) c/ 



9=0J ^ 1 *!(t + 2)!(9-*-2)'- 

Kg)! 



T 2 (Cz) :=I 



12 



9=0 



and 



9-2 



2/ + 1 



,f 0i ti(?-*-2) ! *!(* + 2)!^Z(Z + l) 



are derived in formulas 4 to 7 and 9. Using these closed form formulas, Eq. (3.13) is seen to be identically 
true by simple algebra. □ 

Lemma 6. For every integers p,q>0 the sequences of coefficients defined by the generating functions 
given in Theorem 2 satisfy the recursion 



3 ™ - 5pfl kkTi 







P .\ 



k=0 ' 



P(P+1) 



+ Sp,i 



+ ^ («-l) ! 



n-2 



®q.n-2 
(-if 



(-if 



( n - q -2y( q y(n- q y 

(-\f- l n 



(3.14) 



L"k,q v-> k,1-k-2 yi c k,q-k-2 
k + 2 k + 2 k + 2 



(ky(n-k-2y{k + 2) 2 ( n -\)\(n+\)\ 
C 2 c?„ , ,1 

n> 1. 
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Proof. Applying Lemma 5 and computing the generating functions of both sides of the asserted equality 
we find that the claim of the lemma is equivalent to 

A ^ to A P-M« 



(3.15) 1 -A p , g (z)=8 p ,o z{l :] 
where explicit expressions for 



Ai,,(z)-0 P , 



P(p + i) 
S 2 (z) , S l (z) 



2(1 -Z) 

{-zf- q 



(G(z) - o(z)) + + S 3 ( z ) + U 1 (z) + z + 1 



are derived in formulae 3 and 8. For p = and p > 1 Eq. (3.15) follows immediately from the defining 
equations (2.10b) and (2.10c). For p = 1 the claim follows from combining Lemmas 2 and 3. □ 

Lemma 7. For every integers p,q>0 the sequences of coefficients defined by the generating functions 
given in Theorem 2 satisfy the recursion 



(I 



k,q 



* (-I)"'' , ft 

^ ( n -l)> +Uc '-"- 2 („- <? -2)!( (? )!(»- (? -l)2 



(3.16) 



+ E 

/t=0 



k+1 







P ,i 



h n 

°P = hq 
P(P + 1) 



'q,n—p—l / \| 



) ! (?) ! (' 

(-l)»(2p + l) 



(p)i( B -p- 1)1^+1) 



n> 1. 



Proof. The same as before. We show the equality for every n by showing the equality of the generating 
functions. We obtain that the claim is equivalent to 



=8, 



s ! W , S 2 (z) 



2(1 -z) 



Ai 4 (z) + 5 ( 



!p,o- 



(-z)'' +1 



S 4 (z) + U 2 (z)+z 



(3.17) 



S 2 (z) 
z(l-z] 

where the functions 

S 4 (z):= t 



[G(z)-a(z)]+H(z) 



(-*)"" 



(<z) ! 

B P -!,,(z) (-z)^+'(2 P +l) 

+ (p)\ q y- P {p+i) 

,k+2 it 



and U 2 (z) := £ 



2/ + 1 



„ =(?+ 2(»-<7-2) ! (»-<7-l) 2 ~ v *"' £i(*+l)! 2 fe/(Z + l) 

are evaluated in formulas 10 and 11. For p = Eq. (3. 17) follows from the observation that 

H(z)=z[S 4 (z) + U 2 (z)+z + H(z)]. 
Next we observe that Eq. (2.1 lb) implies 

zB p ^{z) _{-zy + " +2 £ 2k + l {-z) p+q+2P r 2k+l 



B M (z) + 



(-z) p+,?+2 2p + l 



and thus Eq. (3.17) also holds for p > 1. Finally, for = 1 we need to show that 



Bi,,(z) 



2(1 -z) 



Ai,,(z)- 



Bo,,(z) 3(-z 



4( q y 



which, after using the defining equations (2.10) and (2.1 1) several times, amounts to showing that 



S 2 (z) 



1 



.2(l-z) + zJ GW - 
which is exactly what Lemma 2 asserts. 



S'(z) , S 2 (z) | 1 
z 1 — z z 



. . H(z) 3z 2 n 
a(z) + -^ + — =0, 
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Lemma 8. For every integers p,q>0 the sequences of coefficients defined by the generating functions 
given in Theorem 2 satisfy the recursion 



„«+l 

-p,q 



-Spa E 



k=Q 



n" r" 



«-"- 2 \n- q -2y( q y(n- q -iy 



(-i) n 



'n-2 



(-1)" 



- {k)\n-k-2)\k+\) 2 n\ 2 



(3.18) 



n h" I- 1 h" 9 _1 c n 

Y *>9 V" k,q-k-\ yi C k,q-k-\ 

L k4-\ '—' kA- 1 £-1-1 



n> 1. 



Proof. The proof goes along the same lines as the previous ones. Equating the generating functions of both 
sides and using Lemma 2 we need to show that 



C P , g (z) 
z 



=5, 



3 P,0 



z 2(1 -z. 
'(-*)* / S 2 (z) 



4 ^ 

H k=i 



(qf \z(i-z) 

2(- z )« 



(G(z) - a(z)) +H(z) + S 4 (z) + U 2 (z) + 



9 + 1 



q-l 



i-z) 



q+1 



k-iy(k+i)\i ^(ky-{q-k-\yik+\f 



(-z) 



q+1 



k=0 

21 + 1 



(3.19) 



9-1 

"£(9-fc-l) ! (^+l)! 2 ^ 2 '(' + l) 

(_ z )9+l \1 



7 g+l 9-1 



A = l 



2(-l) 



A-l 



(*+l) 



©5,1 Z* +1 4 



9-1 



(9-l) ! 



For p > 1 Eq. (3.19) is immediately clear from the defining equations (2.12b) and (2.12c). For /? = we 
show that the ^-generating functions coincide. Doing this we obtain after some algebra that Eq. (3.18) is 
equivalent to 

S 2 (z)" 



= 
(3.20) 

with the functions 



Z-l 



H(z)+S 4 (z)+z + U 2 (z)-Cz 2 



T'(Cz) 



1 



Cz 2 



5z, 



D(Cz) + ^E 3 (Cz)+z[E 4 (Cz)+T 2 (Cz) 



9-1 



2(-fc)« 



e n /e, (?-^-i) ! (^+i)! 2 ' 



9=0 *=1 

5-1 



z\tz) EE 



9=04=0 



e w ! (?-£-i) ! (£+i) 2 



and 



2/ + 1 



? =o*=i(?-^-l) ! (^+l)! 2 S^ + l) 

given in formulae 12 to 14. Since all functions occurring in Eq. (3.20) are explicitly known the result 
follows from basic algebra. □ 

With this work done we can now prove our main theorem. 

Proof of Theorem 2. The Chapman-Kolmogorov equation implies the recursion 

(3.21) K"{r' ,r)= f ' K{r , ,s)K"- 1 (s,r)ds, r/eR, n>\, 

where K is the one-step transition kernel of A given in Eq. (3.1) From this we can first prove the asserted 
symmetry K" (/, r) = K"(—r', — r) by induction on n. For n = 1 this is clearly true so assuming that it holds 
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for some n > we conclude that K n+l (r 1 ,r) is equal to 

/ K(r',s)K n (s,r)ds = f K(-r',-s)K n (-s,-r)ds 
Jr Jr 

= f K(-r',s)K"{s,-r)ds = K n+l (-r' 1 -r). 



In the next step we prove Eq. (2.9), also by induction on n. For n = 1 the claim is true by Lemma 4. We 
now assume that Eq. (2.9) holds for some n > 1. It then follows that for r > 0: 

K n+1 {r' 7 r) = / K(r',s)K n (s,r)ds 



K(r',s)K"{s,r)ds + / K(r' ,s)K n {s,r)ds + / K(r' ,s)K n (s,r)ds 



o 

n ~0 (iyi-1 -(n+l)r r r 

P ^o J -°° J o 

n-l C_i\n p -(n-p)r , r 

+ E TYT W / 

+ E ^,^" (g+2)r [ r K(r>,s)e-r s ds+ ^ ^ ' f~ K^^^ds 
n-l (_-i\n p -(n-p)r r r <*> 

+ E fvr — 



E c P, q < ri9+2)r [ K(r', S )e-P s d S . 



p,q=0 

The five types of integrals occurring in this expression are easily evaluated to give for /? > 0: 

1 p+l (p+l) (p+2) r - U 



y(r',{)c J ds = <( 1 + r' < r' < r 



— r' > 

P+2 r > u 

r'<0 



r r > r 

P+2 + (p+2) 2 

Jo [ ' j I (p+Tj 2 (p+TH^+2p (p+TKp+2) (p+2) 2 p+^~ U<r<r- 

e -(p+l)r r l_ e -(p+l)r 



P+l ^ (p+1) 2 
.y-(?+2)r 



r' > r 



[ r K(r J ,s)e- ps ds: 
Jo 

J" K(r , ,s)e- ps ds: 



r 1 <0 



p+2 

i _ /-(p+ 2 )' _ gZ (p+iy n , 

p+l p+2 (p+l)(p+2) V<r ^ r ' 

l_ e -(p+l)r , 

p+l r >r 

r < r 



p+2 

-(p+l)r , 

r > r 



p+l (p+l)(p+2) 

This implies that for r 1 < the function K" +l is given by 

(3.22) K" +l (r',r)= £ ^+V / -(« +2 > 

p, 9 =0 
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where 



a" +1 



k~Q k+1 P(P + 1 ) 

+^>q.n-2 



''"-'(n-l)' *■"(«-!)! 



n-2 



(*-«-2)'( 9 )>(n-*)2 
(-1)" 



(3.23) 



' 5 ^{k)\n-k-2)\ P + 2f- 

(-1)" 

t Q k + 2 L k + 2 L k + 2 «•"(«- l)"(n+l) 



-l 



By Lemma 6, a" p + q is equal to a*+ . Similarly, for < r 1 < r the function K" + takes the form 

n n-1 (_\ \n+l J „-p(r' -r)-(n+l)r n+ln+1 

(3.24) K»+\r',r) = £ £ 1 lj , r g _ _ + £ £ &»+V^+ 2 )', 

9 =0 p=0 W ^" ^ p=0o=0 



where 



g«+i _<j 



(-ir 1 (-i> 



i-l n-2 



(-1)" 



(«-l) ! («+l) it to(fe) ! («-fe-2) ! (fe + 2) 2 



c/-2 un q-2 Ji 

^ k,q-k-2 _|_ ^ 2 

fc=0 ^ + 2 ife=0 fc + 2 



and 



(3.25) 



jfe=0 



9 M - "- l L t x 2 U/, ' 1 ^"- p - 1 (^-l)'(n-^-l)! / ,2( p + l)2 



(-l) "(2p+l ) 

(-1)" 



\2 ' i- 



{n-q-lYiqyin-q-lf to k +\ 

Lemma 5 implies that /3^ +1 = Sq^ n V?| and by Lemma 7, is equal to fe^ 1 . Finally, for r' > r the 
function becomes 



(3.26) tf" +1 (r',r) 
where 



(-1)" 



—nr' — r n— 1 / ^\n+l^g—p{r'—r)— (n+\)r «+!«+! 



(I.)' 



p=0 



(p) I (»- J p-l) 1 



p=0o=0 



-pr'-(a+2)r 



■>p,0 



®o,n-2 



(-1)" 



(n-«-2)'(9)'(»-9-l) 2 



km*-k-w+i) 2 k^ 1 k k+i 



(3.27) 



~&q.n- 



n-2 

(-i) 



<:=0 



n-l o-l r " 



"k,q-k-\ 

k k+i 



n 



|2 



<r=0 



In Lemma 8 it was shown that c^ 1 equals c£+ . Combining Eqs. (3.22), (3.24) and (3.26) proves the 
theorem because it follows that for r > the values (/, r) are given by 



i-ip,q=0"'p l q 

(-l)"/-("+ 2 )'- „„-! (-l)"+'/ e -P( r , -0-("+l)-- v „+l v „+l 



r' <0 



-P=0 (p)'(n-p-l) 1 



+ r p t ^=O^V^ /Jr ~ (,+2) ' < r' < r . 



(-ire-" 7 -' , r »-l (-l)"+ 1 /- e -P('- , - r )-'"+ 1 > v „ + l r „+i , + i p /-( 9+2 ) r , 
^ («) ! Lp=0 ( P y (n-p-iy ^ Lp=0Lq=0 c p,q e ' 



> r 
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4. Discussion 

The way in which Theorem 2 is proved gives little insight into how one arrives at the expressions (2. 10) 
to (2.12) for the generating functions A pq , B pq and C Pt q in the first place. It appears to be indicated to 
briefly comment on how we derived these formulas. The first step was to compute the kernels K" for 
low values of n from the Chapman-Kolmogorov equation (3.21) and to observe that they have the form 
asserted in Theorem 2. In the next step we guessed the expression for the part of K" (V, r) not involving the 
coefficients a" p q , b" p q and c" p q so that the problem was reduced to solving the recurrence equations (3.23), 
(3.25) and (3.27). Assuming the validity of Lemma 5 it turns out that the first two of these recurrence 
equation can relatively easily be solved first for G(z), which is — up to the factor {—z) q /{q)' — the 
generating function of (a" + b" )„> i , and then also for A Piq and then B pq . The third recursion for (c n p q ) 
was simplified by the empirical observation that 



q-1 un q-\ _n 

v-< a k,q-k-\ V> c k.q-k-\ _ ~ 

L ~T+] L 1 + i _ q ' n+l 

k=0 K+L k=0 K+i 



i\* ^(uy-iq-k-iyik+iy d * 



for some real numbers d q , q>0, and then solved for C p . q . The educated guesses made in the course of this 
derivation are justified ex posteriori by the proofs presented in this paper. 

Our original motivation was to derive an explicit expression for the asymptotic variance as in Eq. (2.7). 
For this purpose, knowledge of the generating function of the coefficients of the n-step transition kernel, 
as opposed to knowledge of the coefficients themselves, is sufficient. In order to evaluate the infinite sum 
appearing in Eq. (2.7) one is primarily interested in sums of the form Y,n>l a P q which is equal to A pq {\), 
provided this number is finite. Carrying out the computations, however, turns out to be quite subtle and the 
results will be reported elsewhere. 

It is a natural question if the results presented in this paper can be extended to the first-passage percola- 
tion problem on N x {0, 1 , . . . , k}, k > 2. Conceptually, our approach carries over to this setting only if one 
considers semi-directed percolation in which the horizontal edges may be traversed in only one direction; 
the combinatorics involved in computing the one-step transition kernel of the Markov chain A as well as 
the explicit iteration of the Chapman-Kolmogorov equation (3.21), however, soon become unmanageable 
for larger values of k. For the undirected first-passage percolation problem there is the possibility that the 
shortest path {(0,0) = po,pi,. . .,pn-i,Pn = («,0)}, pi — (xi,yi), between (0,0) and (n,0) backtraces by 
which we mean that there exist indices < i < j <N such that Xj < x, . The possible occurrence of such 
configurations prevents an extension of our recursive method to broader graphs in the undirected setting. 
One might also wonder if similar results can be obtained for more general class of edge-weight distributions 
P. It is easy to see that the Markov property of A does not depend on the choice of P and an analysis of 
our proofs shows that the validity of the central limit theorem as well as expression (2.7) for the asymptotic 
variance is not affected by choosing a different edge-weight distribution either, provided one can prove that 
the stationary distribution ft and the one-step kernel K satisfy the moment and mixing conditions used in 
the proof of Theorem 1 . It is however, very difficult, to evaluate the formula for the n-step transition kernel 
explicitly, if P is not the exponential distribution, although our approach via generating functions remains 
likewise applicable. 
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